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INTRODUCTION 

In  the  last  two  decades,  a  tremendous  amount  of  work  has 
been  done  to  improve  methods  of  neutron  detection.   The  behavior 
of  neutrons  in  matter  is  quite  different  from  that  of  either 
charged  particles  or  gamma  rays.   Since  the  neutrons  are  un- 
charged, no  Coulomb  forces  come  into  play  with  either  the  orbital 
electrons  or  the  nuclei.   Thus  for  neutrons  to  affect  matter, 
they  must  either  enter  the  nucleus  or  come  sufficiently  close  to 
it  for  the  nuclear  forces  to  act.   There  are  several  known  mechan- 
isms by  which  the  interaction  of  neutrons  with  matter  takes  place 
and  each  of  these  processes  is  the  basis  of  a  potential  method  of 
detection.   For  example,  in  the  neutrons -induced  transmutations 
(n,  oC),    (n,  p),  and  (n,  fission)  reactions,  the  product  par- 
ticles make  possible  the  detection,  while  in  transmutations  which 
result  in  radioactive  product  nuclei  the  subsequent  decay  of  the 
radioactive  nuclei  gives  information  on  the  neutron  flux  which 
induces  radioactivity. 

Among  the  most  extensively  used  methods  is  the  B^  (nt<X) 
reaction.   Many  factors  make  it  highly  satisfactory  for  this  pur- 
pose, namely,  the  high  cross  section  of  the  reaction  and  the 
simple  energy  dependence  over  a  wide  range  of  values.   The  reac- 
tion is  easy  to  detect  even  In  the  presence  of  comparatively 
large  gamma  fluxes  because  of  the  high  specific  ionization  and 
the  large  energy  of  the  charged  particles  which  are  released. 
Solid  boron  as  well  as  gaseous  compounds  can  be  used,   boron  tri- 
fluoride  (BF3),  being  a  comparatively  stable  compound,  is  used 
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frequently.  Proportional  counters  filled  with  BF3  gas  have 
played  an  important  part  in  the  development  of  nuclear  technol- 
ogy. They  have  been  widely  applied  in  the  detection  of  thermal 
and  epithermal  neutrons  in  a  variety  of  forms  that  range  from 
multiple  arrays  of  large  counters  providing  high  efficiency  de- 
tection, for  example,  in  time  of  flight  spectrometers,  to  small 
counters  which  permit  determination  of  flux  distribution  in  re- 
actor core  assemblies. 

The  purpose  of  this  investigation  was  as  follows* 

1.  To  get  a  better  understanding  of  the  different 
problems  encountered  in  the  construction  and  preparation  of 
boron  trifluoride  counters. 

2.  To  prepare  the  best  BP3  counter  possible  by  using 
an  inexpensive  vacuum  system. 

3.  To  study  the  characteristics  of  the  BP3  counter 
obtained. 

4.  To  use  the  prepared  counters  to  determine  certain 
subcritical  reactor  parameters. 

PREPARATION    OF  THE   COUNTERS 

Construction  of  the  Vacuum  System 

Since  the  production  and  purification  of  boron  trifluoride 
gas  had  to  be  done  under  vacuum,  a  special  design  of  vacuum  sys- 
tem was  necessary.   The  vacuum  system  had  essentially  to  perform 
two  functions: 


1.  Evacuate  the  whole  system,  including  the  counters. 

2.  Purify  the  gas  under  vacuum. 

The  photographs  of  the  vacuum  system  are  shown  in  Pigs.  1 
and  2.   A  Duco-3eal  pump,  B,  with  an  exhaust  vent,  I,  was  used 
to  pump  the  pressure  down  to  five  microns.   Mechanical  vacuum 
pumps  remove  gases  from  the  system  by  compressing  them  into  a 
smaller  volume,  thus  increasing  their  pressure  until  they  force 
an  oil- sealed  valve  open  and  then  escape  into  the  atmosphere. 
This  compression  works  very  well  with  the  permanent  gases  such 
as  air  and  hydrogen,  but  condensable  vapors  like  water  vapor  and 
solvent  vapors  may  cause  a  great  deal  of  trouble.   These  vapors 
raise  the  vapor  pressure  of  the  oil,  and  thus  limit  the  vacuum 
attainable  by  the  pumps.   The  problem  of  obtaining  high  vacuum 
in  presence  of  the  condensable  gases  was  solved  by  introducing 
an  exhaust  vent  in  the  pump.   A  valving  arrangement  is  provided 
which  permits  atmospheric  air  to  enter  the  exhaust  cycle  of  the 
pump.   Because  of  the  addition  of  this  extra  volume  of  air,  the 
vapors  coming  from  the  system  can  be  compressed  only  slightly  to 
force  the  exhaust  vent  in.   Thus  since  the  vapors  are  not  com- 
pressed to  the  critical  point,  they  do  not  condense  but  pass  in- 
stead in  vapor  form  out  of  the  pump.   Even  if  some  condensable 
material,  such  as  water,  should  get  into  the  oil  in  liquid  form, 
the  vented  exhaust  will  clear  the  oil  in  a  matter  of  seconds. 
This  exhaust  vent  has  proven  to  be  effective  even  in  the  presence 
of  BP3  gas. 

A  Consolidated  Electrodynamics  glass  oil  diffusion  pump,  A, 
Pig.  1,  was  used  to  reduce  the  pressure  to  the  range  10   to  10 


Fig.  1.   Photograph  of  vacuum  pumps  and  associated  control  panels 
for  the  ionization  and  Pirani  gauges.   A,  diffusion  pump;  B, 
mechanical  pump;  C,  Pirani  gauge;  D,  control  panel  RG-21; 
E,  control  panel  for  GP-110  (Pirani  gauge);  P,  potentiometer; 
G,  variac;  H,  VTVM  Heathkit;  I,  exhaust  vent. 
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mm  Hg.   It  consisted  of  three  independent  boilers,  each  contain- 
ing a  heater  immersed  directly  in  the  operating  liquid. 

This  pump  was  mounted  on  adjustable  clamps  and  connected  to 
the  system.   The  whole  system  was  then  tested  by  a  high-frequency 
Tesla  coil.  When  all  the  leaks  were  sealed,  the  diffusion  pump 
was  filled  with  Dow  Corning  diffusion  pump  oil  type  703  to  ap- 
proximately five  mm  above  the  heaters.   After  obtaining  the  back- 
ing pressure  from  the  mechanical  pump,  the  heaters  were  connected 
to  the  a-c  mains  through  a  variac  transformer,  G.   Prolonged  de- 
gassing of  the  oil  was  necessary  to  obtain  the  lowest  pressure. 
In  general,  twelve  to  fourteen  hours  were  necessary  to  obtain 
pressure  of  the  order  of  10'5  mm  Hg.   The  right  voltage  for  the 
heater  was  selected  by  operating  the  diffusion  pump  at  different 
voltages  and  checking  the  pressure  in  the  system  with  the  ioniza- 
tion gauge  control,  D.   The  lowest  pressure,  5  x  10~6  mm  Hg,  was 
obtained  at  70  volts.  This  was  then  selected  as  the  operating 
voltage.  A  general  precaution  is  to  avoid  putting  the  diffusion 
pump  to  work  before  the  pressure  is  reduced  down  to  10  microns. 
For  this  reason  a  Pirani  gauge,  C,  type  GP-110,  obtained  from 
Electrodynamics  Corporation,  was  Installed  between  the  mechanical 
pump,  B,  and  the  oil  diffusion  pump,  A.   Pressure  in  the  range  of 
0  to  2000  microns  could  be  measured  by  this  gauge. 

For  measuring  pressures  below  one  micron,  Veco  ion  gauge 
type  RG-75  together  with  a  control  panel,  D,  type  RG-21,  was 
used.   The  gauge  contains  a  thoria-coated  indium  filament  which 
does  not  oxidize  if  the  filament  accidentally  comes  in  contact 
with  air.   The  RG-75  was  connected  directly  to  the  manifold,  F, 


as  shown  in  Pig.  2.   In  a  later  design,  the  position  of  the  gauge 
was  changed  as  seen  in  Fig.  9,  because  it  was  found  that  a  thin 
white  film  settled  on  its  walls  and  changed  its  characteristics. 
Heating  the  tube  walls  to  about  250°  C  helped  in  removing  that 
film.   The  tube  grid  had  to  be  outgassed  at  a  pressure  less  than 
10~4  ram  Hg,  because  outgassing  at  higher  pressures  tends  to 
blacken  the  walls  due  to  sputtering  of  the  tungsten  in  the  grid. 

Two  Alpert  valves  V^  and  V^,  Pig.  2,  were  used  to  substitute 
for  stopcocks  between  the  pumps  and  the  manifold  and  to  discon- 
nect the  counters  from  the  rest  of  the  system.   This  type  of 
vacuum  valve  was  developed  at  the  Westinghouse  Laboratories.   It 
essentially  consists  of  a  separate  valve  assembly  which  becomes 
a  part  of  the  vacuum  system,  and  a  driving  mechanism  which  opens 
and  closes  the  valve.   Aside  from  the  vacuum  leads,  the  valve 
assembly  contains  only  three  mechanical  parts:   A  copper  cup, 
1.75-inch  diameter,  a  flexible  Kovar  diaphragm  which  allows  a 
movement  (up  and  down)  of  about  0.1  inch,  and  a  Kovar  nose,  with 
a  highly  polished  45-degree  conical  surface,  seated  in  the 
copper  cup. 

All  of  the  joints  in  this  valve  were  brazed  (eutectic  silver- 
copper  soldered)  simultaneously  in  a  hydrogen  furnace  without 
the  use  of  flux.   The  driving  mechanism  consisted  of  an  outside 
screw  and  housing  (stainless  steel),  a  simple  screw  (silicon 
bronze),  and  a  drive  screw  (stainless  steel)  which  was  bolted  to 
the  valve  nose  before  the  rest  of  the  driver  was  installed.   The 
external  portion  of  the  nose  was  machined  to  form  a  flat  disk 
one  inch  in  diameter.   The  simple  screw  (32  threads  per  Inch) 


IONIZATION    GAUGE 
COUNTERS 


3     ft  A- 


Fig.  2.   Photograph  of  the  filling  system.   A,  50-ml  flask; 
T3>  T4^  T5,  Tg,  traps;  V-^,    V2,  Alpert  valves;  M,  manom- 
eter; P,  manifold. 
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exerts  its  thrust  against  this  flat  disk  which  forms  the  thrust 
bearing.   During  the  bake-out  periods,  the  driving  mechanism  was 
removed  completely. 

A  trap,  Tg,  Pig.  2,  similar  to  the  one  used  by  Alport  (2), 
was  placed  between  the  high  vacuum  side  and  the  diffusion  pump 
to  prevent  the  backs treaming  of  oil  vapor  in  the  high  vacuum 
side  of  the  system.   The  trap  consisted  of  a  standard  reentrant 
glass  trap  into  which  corrugated  copper  was  inserted  so  that  the 
trap  contained  a  large  number  of  fine  copper  straws.   This  type 
of  trap  was  designed  originally  to  solve  the  usual  problem  en- 
countered with  a  refrigerated  reentrant  trap,  that  is,  as  the 
liquid  refrigerant  evaporates,  the  cold  level  drops  and  the  gas 
that  had  previously  been  absorbed  is  reevaporated  into  the  system. 
In  this  design,  it  was  intended  to  utilize  the  high  thermal  con- 
ductivity of  the  metallic  insert  so  that  the  entire  trap  surface 
may  remain  cold  despite  large  fluctuations  in  the  liquid  level. 
The  trap  also  affords  an  extremely  high  ratio  of  length  to  diam- 
eter without  greatly  sacrificing  the  speed  of  the  system.   It 
was  found  later  on,  (3),  that  the  trap  worked  efficiently  even 
when  the  liquid  level  dropped  to  zero.   Since  then,  such  traps 
have  been  used  in  the  nonreentrant  portions  of  vacuum  systems 
without  refrigerants. 

A  mercury  manometer,  M,  was  used  to  measure  the  BF3  gas 
pressure  in  the  system  during  the  process  of  filling  the  counters. 
A  layer  of  703  Dow  Corning  oil  was  placed  on  top  of  the  mercury 
in  order  to  avoid  the  Hg  pressure  which  amounts  to  1.3  x  10"4  mm 
at  room  temperature.   About  four  to  six  hours  were  necessary  to 


pump  down  all  the  entrapped  gases  from  the  oil.   Long  contact  of 
BP3  gas  with  the  oil  seemed  not  to  affect  it  to  any  appreciable 
extent. 

A  water  pump  was  connected  to  the  system  by  high  vacuum 
rubber  tubing  through  a  phosphorous  pentoxide  trap  to  remove 
the  BF3  gas  from  the  system  when  needed. 

As  the  work  proceeded,  the  BP3  gas  was  found  to  react  se- 
verely with  the  stopcock  grease  introducing  impurities  in  the  gas 
and  spoiling  the  counter  characteristics.   It  was  decided  to  mod- 
ify the  system  by  cutting  down  the  number  of  glass  stopcocks  to 
one.   The  modified  system  is  illustrated  in  Pig.  9. 

A  Cu-Constantan  thermocouple  was  also  placed  inside  trap  T3. 
This  thermocouple  was  used  to  measure  the  temperature  of  the  gas 
evolved.   The  Cu-Constantan  wires  were  sealed  through  an  Elec- 
tric Industries  strain-free  seal  type  D-60  TT-SX,  which  was  soft 
soldered  to  the  Kovar  glass  seal,  as  shown  in  Pig.  9.   The  glass 
side  of  this  seal  was  connected  to  trap  3.   The  end  of  the  ther- 
mocouple was  placed  about  one  inch  above  the  bottom  of  the  trap. 
A  precision  potentiometer  type  K3,  Leeds  and  Northrup,  was  used 
to  measure  the  potential  difference  across  the  wires.   Tempera- 
tures were  determined  by  referring  to  the  appropriate  calibrated 
charts. 

Construction  of  the  Counters 

Two  fundamental  principles  should  be  considered  in  the  con- 
struction of  the  counters: 
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1.  The  construction  and  preparation  should  be  such  as 
to  permit  the  mechanism  of  the  counter  discharge  to  operate  with 
a  minimum  of  disturbing  effects  from  unwanted  agencies. 

2.  The  construction  procedure  should  insure  that  the 
desired  conditions  remain  as  permanently  as  possible. 

The  first  principle  requires  that  the  counter  cylinder  should 
be  made  of  a  substance  with  high  work  function,  and  alkali  metals 
should  be  especially  avoided.  This  is  important  because  the  sup- 
ply of  secondary  electrons,  which  causes  the  discharge  to  con- 
tinue and  thus  defeat  the  processes  quenching  the  discharge,  is 
dependent  on  the  work  function.   The  number  of  secondary  elec- 
trons can  be  greatly  reduced  by  proper  choice  of  cathode  surface, 
and  the  counter  will  accordingly  have  better  operating  charac- 
teristics. 

The  second  principle  requires  cleaning  the  counter  thoroughly 
before  sealing  to  the  vacuum  system.   Any  absorbed  film  of  ma- 
terial will  change  the  chemical  composition  of  the  gas  in  the 
counter,  thus  changing  the  properties  of  the  counter  accordingly. 
This  is  especially  true  with  the  counters  filled  with  boron  tri- 
fluorlde  gas. 

Copper  has  been  found  to  be  a  good  material  for  the  counter 
cylinders.   Brass  could  also  be  used  and  is  frequently  available 
in  suitable  shapes.   Both  these  materials  were  used  in  this  in- 
vestigation and  the  counters  built  are  illustrated  in  Pigs.  3 
and  4. 

The  outer  shell  of  counter  (A),  which  served  also  as  the  cath- 
ode, was  made  of  a  seamless  brass  tube  seven  inches  long,  one 
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inch  outside  diameter,  with  a  wall  thickness  of  0.03  inch.   The 
anode,  1,  was  a  3-mil  diameter  tungsten  wire.   Tungsten  was 
selected  because  it  could  be  readily  sealed  to  glass  and  heated 
successfully  while  the  counter  was  being  pumped.   A  piece  of  uni- 
form, smooth,  and  clean  tungsten  wire  was  selected  to  avoid  the 
local  ionization  which  produces  spurious  counts,  and  to  avoid 
the  spread  in  gas  multiplication  and  the  spread  in  pulse  height 
caused  by  the  small  eccentricities  in  the  wire  (21).   The  end  of 
the  anode  wire,  1,  was  spotwelded  at  point  2,  to  the  anode  lead 
3,  which  was  sealed  to  glass  at  point  4.   Point  2,  where  the  wire 
was  welded  to  the  supports,  was  shielded  with  a  glass  sleeve,  5, 
in  order  to  avoid  the  high  local  fields. 

Details  of  the  end  cap,  6,  are  shown  in  Pig.  3(B).   It  was 
made  of  a  solid  piece  of  copper  cylinder  one  inch  in  diameter  and 
5/8  inch  long  through  which  a  cylindrical  hole  1/2  inch  in  diam- 
eter and  9/16  inch  long  was  drilled.   A  l/2-inch  diameter  quartz 
tube,  7,  5/8  inch  long,  with  a  closed  end,  was  placed  in  the 
drilled  hole.   A  small  hole  was  blown  in  the  closed  end  by  a 
fine  oxygen -hydrogen  flame  before  the  quartz  tube  was  placed  in 
the  drilled  cylindrical  hole.   The  anode  wire,  1,  was  introduced 
through  the  hole,  through  the  Nichrome  spring,  8,  and  through  a 
7/16-inch  diameter  molybdenum  disk,  9,  in  which  a  fine  hole  had 
been  drilled.   Pressure  was  applied  on  this  disk  and  a  small  bead 
of  pyrex  glass,  10,  was  sealed  to  the  wire  to  keep  it  under  ten- 
sion.  A  l/4-inch  copper  disk,  11,  was  silver  soldered  as  shown 
in  Pig.  3(B)  to  close  the  end.  A  glass-copper  seal,  12,  Pig.  3 
(A),  was  used  as  an  insulator  for  the  high-voltage  electrode  and 
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also  for  connecting  the  counters  to  the  vacuum  system.   Cylinder 
13  is  a  piece  of  brass  tubing  one  inch  in  diameter  and  1/8  Inch 
thick  which  was  silver-soldered  to  the  outer  shell  of  the 
counter.   A  copper  disk,  14,  3/32  inch  thick,  in  which  a  5/16- 
inch  hole  was  drilled,  was  seated  in  the  cylinder,  13,  and  silver- 
soldered  to  It.   A  glass-copper  seal,  12,  was  fitted  by  force 
into  this  hole  and  was  silver -soldered  to  the  disk,  14.  Details 
of  this  part  of  the  counter  are  shown  in  Fig.  3(C). 

Counter  (B),  Pig.  4,  was  constructed  from  a  special  1-inch 
diameter  glass-copper  seal,  1,  obtained  from  Larson  Electronics 
Corporation,  California.   The  copper  part  of  this  seal  was  7 
inches  long.  The  end  cap,  6,  was  the  same  as  the  one  used  in 
counter  (A).   This  counter  was  easier  in  construction  and  proved 
to  be  better  in  quality,  and  was  therefore  preferred  in  most  of 
the  work. 

Before  connecting  the  counters  to  the  vacuum  system,  they 
were  cleaned  with  dilute  nitric  acid,  washed  thoroughly  with 
distilled  water,  and  dried  in  an  oven  at  200°  C,  thus  removing 
most  of  the  absorbed  impurities.   After  connecting  the  counters 
to  the  vacuum  system  they  were  outgassed  at  250°  C  for  12  hours 
and  then  at  400°  C  for  two  more  hours.   The  pressure  during  the 
bake -out  period  was  kept  below  10"4  mm  Hg.   Baking  out  of  the 
counter  envelopes  proved  to  have  a  great  influence  on  the  counter 
behavior  later  on. 
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Effect  of  Impurities  in  the  Boron  Trifluoride 
Gas  on  the  Counter  Performance 


On  account  of  the  unsatisfactory  behavior  of  commercial 
boron  trifluoride  gas  as  a  direct  filling  for  pulse  counters  of 
high  sensitivity,  some  preliminary  work  was  carried  out  at  the 
Chalk  River  Laboratories,  in  Canada,  on  the  purification  of  boron 
trifluoride.   It  was  found  that  the  impurities  most  likely  to 
interfere  with  the  performance  of  boron  trifluoride  counters  are 
oxygen,  halogens,  ammonia,  water,  nitrous  and  nitric  oxides,  sul- 
phur dioxide,  hydrogen  fluoride  and  chloride,  silicon  tetraflu- 
oride,  and  mercury  vapors.   Their  influence  has  been  attributed 
to  electron  capture  and,  in  some  cases,  to  dissociation  of  the 
molecules  (22). 

The  electron  attachment  probability,  h,  for  some  of  the 
gases  mentioned  above  are  presented  in  Table  1,  (16).  It  is 
defined  as  the  probability  that  an  electron  will  attach  to  the 

molecule  with  which  it  collides.   An  electron  will  make,  on  the 

5 
average,  around  10  collisions  In  crossing  a  counter,  the  mean 

free  path  depending  on  the  pressure.   Considering  this,  the 
probability  that  the  electron  should  attach  to  a  neutral  mole- 
cule should  be  less  than  10~5.   Reported  values  for  h  are  low 
for  argon,  carbon  monoxide,  hydrogen,  nitrogen,  and  ammonia  (16). 
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Table  1.  Electron  attachment  probability,  h  (16). 


Gas     : 

h        : 

Gas 

:         h 

CO 

6.25  x  10-7 

H20 

0.25  x  10"5 

NH3 

1.01  x  10~6 

ci2 

0.5  x  lO"3 

C2H6 

4.0  x  10"5 

SiP4 

1.48  x  10~5 

N20 

1.63  x  10~4 

so2 

10-4 

Air 

0.5  x  lO"5 

sp6 

0.936 

o2 

0.25  x  10-3 

h  -  probability  that  an  electron  will  attach  to  the 
molecule  with  which  it  collides. 

The  h  for  air  is  0.5  x  10~5,  which  Is  close  to  the  critical 
value  of  10"5,  mentioned  above.   The  figures  for  oxygen,  water, 
and  chlorine  are  all  greater  than  this  critical  value,  and  there- 
fore must  be  avoided. 

Boron  trifluoride  generated  inside  a  glass  vacuum  system 
usually  contains  silicon  tetraf luoride  (6,  14),  an  electronega- 
tive gas  with  attachment  probability  of  1.485  x  10"5.   Sulphur 
dioxide  (7)  and  sulphur  hexaf luoride  (13)  have  also  been  found 
in  BPg.   Both  are  highly  electronegative  gases  with  attachment 
probabilities  of  10"4  and  0.936,  respectively. 

The  effects  which  silicon  tetrafluoride,  sulphur  dioxide, 
and  sulphur  hexaf luoride  have  on  the  plateau  and  the  pulse  size 
distribution  of  a  proportional  counter  were  investigated  by 
Aponte  and  Korff  (4).   They  obtained  the  attachment  cross  sec- 
tions for  these  gases  as  follows: 

The  number  of  electrons  at  a  distance  dr  Is  expressed  as 
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dn  =  io(   -  7"!  )  ndr 
where  (A   is  the  first  Townsend  coefficient,  and  77  is  the  number 
of  attachments  per  cm  advanced  in  the  direction  of  the  field. 
Integrating,  the  number  of  electrons  at  a  distance  r  can  be 
written  as 

n  =  n0  exp  f(OC  -  Y\  )  dr  (2) 

The  number  of  attachments  per  cm  advanced  can  be  written  as 

=  (c/A)(l/v)(p/P)h  (3) 

where  c  is  the  average  agitational  velocity  of  the  electron,  \ 
is  its  collision  mean  free  path,  (c"/A)  is  the  frequency  of  the 
collision,  and  v  is  the  drift  velocity  along  the  field  of  the 
counter,  making  (l/v)  the  time  needed  to  cover  one  cm  along  the 
field.   Therefore  (c/A)(l/v)  is  the  number  of  collisions  the 
electrons  make  while  covering  a  distance  of  one  cm  along  the 
field.   Since  only  a  fraction  of  these  collisions  are  with  the 
electronegative  molecules,  the  number  of  collisions  with  the  im- 
purity can  be  obtained  by  multiplying  the  above  expression  by 
(p/P)»  the  ratio  of  the  pressure  of  the  electronegative  gas  to 
that  of  BF3.   Multiplying  (c/ A) ( l/v) (p/P)  by  h  gives  the  number 
of  attachments  per  cm  along  the  path  of  the  electron. 

Prom  equation  (3),  it  can  be  seen  that  Yl    is  a  function  of 
p,  and  since  c<  is  a  constant  for  the  counter,  by  increasing  the 
pressure  of  the  electronegative  gas,  fewer  electrons  will  reach 
the  wire.   This  will  cause  a  variation  in  the  pulse  size  distri- 
bution of  the  counter,  with  the  consequent  loss  of  the  counting 
rate  when  the  pulse  becomes  smaller  than  the  discrimination  level 
of  the  amplifier.   This  will  affect  the  plateau  of  the  counter. 
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Since  (c/70(l/v)  =  N0,  where  NQ  is  the  total  number  of  colli- 
sions  the  electron  makes  in  going  from  cathode  to  point  r,  then 
equation  (3)  can  be  written  as 

n  =  n0  exp   joCdr  -  h  NQ(p/P)r  (4) 

where  h  was  assumed  as  a  constant  of  the  electronegative  gas. 
If  the  collision  cross  section  is  given  by 

crc   =  1/X  n  (5) 

where  7^  i9  the  mean  free  path  for  collision  of  the  electrons  in 
the  electronegative  gas,  and  N  is  the  number  of  electronegative 
molecules  per  cm  ,  then  the  attachment  cross  section  is  given 

by 

o\  -  crc  h  =  hA »  (6) 

The  effect  of  silicon  tetraf luoride,  sulphur  hexaf luoride, 
and  sulphur  dioxide  on  the  counter  plateau  is  shown  in  Pigs.  5, 
6,  and  7,  respectively  (4).   From  Pig.  5,  it  is  seen  that  a  very 
small  amount  of  silicon  tetraf luoride  is  needed  to  damage  a 
counter.   An  amount  of  silicon  tetraf luoride  equal  to  0.2  per 
cent  reduces  the  length  of  plateau  by  50  volts.   The  attachment 
cross  section  for  this  gas  obtained  by  Aponte  and  Korff  (4)  was 

on    q 

5.12  x  10  *u  cm6.   The  strong  effect  of  sulphur  hexaf luoride  is 
due  to  its  large  attachment  probability.   The  attachment  cross 
section  was  found  (4)  to  be  4.26  x  10" 15  cm2. 

Sulphur  dioxide  also  has  a  high  attachment  probability  of 
10"  .   The  amount  of  this  gas  permitted,  without  the  counter  be- 
ing affected  beyond  the  limits  of  tolerance,  is  given  as  0.01 
per  cent  (4). 

No  direct  measurements  as  to  the  effects  of  these  impurities 
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on  the  counter  performance  were  conducted  during  this  investiga- 
tion, although  various  precautions  necessary  to  avoid  these  im- 
purities were  taken  in  order  to  get  a  counter  with  good  charac- 
teristics.  Stopcock  greases  containing  silicon  or  sulphur  were 
specifically  avoided. 

Preparation  of  Boron  Trifluoride  Gas 

The  nominal  composition  of  commercial  boron  trifluoride  is 
20  atom  per  cent  B10  and  80  atom  per  cent  B11.   The  reaction  B 
(n,  (X)  occurs  only  for  the  B10  isotope.   The  B11  component  per- 
forms no  useful  function.   An  extensive  amount  of  work  on  the 
preparation  of  pure  boron  trifluoride  gas  from  commercial  cylin- 
der gas  was  undertaken  at  A.E.R.E.,  Harwell,  by  Hudswell,  Nairn, 
and  Wilkinson  (14).   Various  precautions  were  necessary  as  boron 
trifluoride  gas  reacts  very  readily  with  the  stopcock  grease  and 
the  silica  in  the  glass.  All  the  stopcocks  were  replaced  by  mer- 
cury cut-off  valves  and  the  glass  container  used  for  storing  BFg 
by  stainless  steel  in  order  to  avoid  any  prolonged  contact  with 
the  glass,  and  polytetrafluoroethylene  was  chosen  as  the  gasket 
material  for  this  vessel.   They  reported  that  rough  fractionation 
of  the  cylinder  gas  by  repeated  condensation  at  -120u  C  and  pro- 
longed pumping  of  the  permanent  gases  from  solid  boron  trifluor- 
ide were  not  effective  in  improving  the  gas  characteristics. 
Hence  they  proceeded  via  the  additives  with  diethyl  ether  and 
prepared  boron  trifluoride  diethyl  ether  complex  by  passing  dried 
boron  trifluoride  gas  through  a  flask  filled  with  analytically 
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pure  ether.   The  product  was  fractionally  distilled  in  an  effi- 
cient column  packed  with  glass  helics  and  the  fraction  boiling 
at  126°  C  was  collected.   They  recommended  fractionation  under 
vacuum  in  order  that  the  product  could  be  stored  for  a  longer 
period  of  time.   Their  next  step  was  the  preparation  of  the 
BF-zCaFg  complex. 

They  dried  CaF2  in  vacuum  at  300°  C  and  added  boron  triflu- 
oride  etherate  to  the  flask.   The  etherate  reacts  with  calcium 
fluoride  at  temperatures  below  100°  C  to  form  CaF2BF3  with  the 
liberation  of  ether.   They  recommended  dry  steel  containers  for 
the  storage  of  this  complex.   The  replacement  of  the  diethyl 
ether  boron  trifluoride  complex  by  the  dimethyl  ether  complex 
was  recommended  by  Lloyd,  Hodges,  and  Weinstock  (17)  because  of 
its  higher  volatility. 

Air-free  boron  trifluoride  gas  was  prepared  from  the  BF^CaFg 
complex  by  Bistiline  (5)  and  later  on  by  Fowler  and  Tuncliffe 
(10).   Their  method  was  adopted  in  this  investigation  with  minor 
changes.   The  detailed  method  used  was  as  follows. 

A  known  amount  of  96  per  cent  enriched  CaF2BF3  complex,  ob- 
tained from  Oak  Ridge  National  Laboratory,  was  kept  in  the  50-ml 
round  bottom  flask  A,  Fig.  8.   After  outgassing  the  counters  as 
described  in  the  section  on  construction  of  the  counters,  the 
system  was  pumped  down  to  less  than  10"**  mm  Hg.   The  Alport  valve 
V^  was  closed,  a  heating  mantle  was  placed  around  flask  A,  and 
the  complex  was  heated  to  100°  C  for  two  hours,  or  until  the 
pressure  was  reduced  to  10  •  mm  Hg,  indicating  the  removal  of 
any  moisture.   Traps  T2  and  T3  were  surrounded  by  Dewar  flasks 
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filled  with  solid  CO2  -  acetone  mixture  and  liquid  air,  re- 
spectively. 

The  complex  was  heated  to  150°  C  and  the  pumping  was  con- 
tinued until  the  pressure  was  less  than  10"°  mm  Hg.   Stopcock  S^ 
was  closed  and  the  temperature  of  the  compound  was  raised  to  be- 
tween 180°  and  200°  C.   The  gas  started  to  evolve.   As  indicated 
by  Hudswell,  Nairn,  and  Wilkinson  (14),  the  gas  evolved  during 
the  first  ten  minutes  included  most  of  the  volatile  impurities. 
For  this  reason  trap  T^  was  cooled  down  by  surrounding  it  with  a 
Dewar  flask  filled  with  liquid  air  and  the  gas  evolved  during 
the  first  ten  minutes  was  collected  in  it.   T^  was  then  sealed 
at  the  constriction.   Alpert  valve  Vg  was  then  closed  with  a 
torque  wrench  and  the  stopcock  S^  was  slowly  opened.   The  gas 
passing  through  trap  T2  was  purified  from  hydrogen  fluoride  (HP), 
which  condensed  at  dry-ice-acetone  temperature,  and  was  collected 
in  trap  T3  at  liquid  air  temperature  in  a  solid  form.   When  enough 
BF3  was  collected  in  T3,  stopcock  S^  was  closed  and  Vg  was  opened. 
The  system  was  pumped  down  to  less  than  10~5  mm  Hg,  V2  was  closed, 
and  the  liquid  air  surrounding  T3  was  removed  to  evaporate  the 
collected  solid  BF-j.   When  the  BF3  pressure,  as  indicated  by 
manometer  M,  acquired  a  predetermined  value,  V]_  was  closed.   The 
remaining  gas  in  the  system  was  removed  by  opening  stopcock  S3 
and  cooling  T4  down  to  liquid  air  temperature. 

The  counters  now  filled  with  BF3  were  heated  for  two  hours 
by  an  infrared  lamp  to  accelerate  the  reaction  of  BF3  gas  with 
the  counter  walls.   This  gas  was  then  removed  from  the  counters 
by  opening  V-j^  and  stopcock  S3  and  collecting  it  in  trap  T4.   The 
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system  was  once  again  pumped  down  to  less  than  10"  5  mm  Hg.   This 
procedure  of  filling,  removing,  and  pumping  was  repeated  three 
times  before  the  final  filling  of  the  counters,  after  which  the 
valve  V^  was  closed  and  the  counters  were  wrapped  by  wet  asbestos 
sheets  and  sealed  at  the  constrictions. 

The  procedure  adopted  for  filling  the  counters  using  the 
modified  vacuum  system  shown  in  Fig.  9  was  as  follows. 

1.  The  counters  were  evacuated  and  outgassed  as  be- 
fore and  the  system  was  checked  for  leaks. 

2.  The  compound  was  heated  to  100°  C  and  the  pumping 
was  continued  until  the  pressure  was  down  to  less  than  10"°  mm 
Hg. 

3.  Alpert  valve  V,  was  closed  and  the  compound  was 
heated  to  between  180°  and  200°  C.   The  gas  evolved  during  the 
first  ten  minutes  was  condensed  in  trap  T^  and  T^  was  sealed  at 
the  constriction. 

4.  Dewar  flasks  containing  solid  COg  -  acetone  mixture 
and  liquid  air  were  placed  around  Tg  and  T^,  respectively. 

5.  Alpert  valve  Vg  and  stopcock  S^  were  closed  and  V-^ 
was  opened.   The  BF3  evolved  was  collected  in  trap  T3. 

6.  When  the  desired  amount  of  gas  was  collected,  V^ 
was  closed  and  the  system  was  opened  to  the  pumps  by  opening 
stopcock  S^. 

7.  The  system  was  pumped  down  until  the  pressure  be- 
came less  than  10"°  mm  Hg.   S^  was  closed  and  the  counters  were 
opened  to  the  system  by  opening  the  valve  V2.   The  liquid  air 
flask  surrounding  trap  T3  was  slowly  removed  and  the  pressure  of 
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the  gas  evolved  was  checked  on  manometer  M. 

8.  Vg  was  closed  and  the  counters  were  heated  by  infra- 
red lamp  for  two  hours.   The  gas  in  the  system  was  collected  in 
trap  T4  which  was  then  sealed  at  the  constriction. 

9.  This  procedure  of  filling,  soaking,  and  pumping  out 
the  BP3  gas  was  repeated  once  more,  before  sealing  the  counters 
at  the  constriction.   Every  time  the  system  was  pumped  down  be- 
fore filling,  the  temperature  of  the  solid  BP3  was  checked  with 
the  Cu-Constantan  thermocouple  located  in  T3.   The  temperature 
was  usually  irregular  for  a  period  of  time  and  then  became  con- 
stant.  This  could  be  attributed  to  the  condensable  impurities 
being  pumped  out.   The  filling  process  was  never  started  until 
this  constant  temperature  was  acquired.   The  temperature  of  the 
evolved  gas,  after  removal  of  the  liquid  air  flask,  was  also 
checked  to  test  the  purity  of  the  gas.   Figures  10  and  11  show 
the  variance  of  temperature  (T)  with  time  (t)  at  two  different 
pressures.   Data  used  in  plotting  the  above  curves  are  listed 

in  Table  2. 

COUNTER  INSTRUMENTATION 

A  block  diagram  indicating  the  instrumentation  generally  re- 
quired for  the  boron  trifluoride  proportional  counter  is  shown 
in  Fig.  12. 

As  the  tube  is  normally  located  at  a  distance  from  the  in- 
strumentation, coaxial  cable  must  be  employed  to  transmit  the 
pulses  from  the  tube  to  the  amplifier.   For  lengths  of  cable  less 
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Table 

2. 

Variation  of 

BF3 

gas 

temperature 

with  time. 

BP3  pressure 

=  12.6  inches 

Hg  : 

BF3  pressure  = 

11.79  inches 

Hg 
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:    T°  C 

• 
• 

t  minutes 

:    T°  C 

0 

-78.4 

0 

-78.3 

5 

-70.0 

10 

-61.0 

10 

-61.4 

15 

-50.5 

15 

-55.0 

20 

-48.0 

20 

-42.3 

25 

-42.5 

25 

-32.6 

30 

-39.0 

30 

-29.0 

35 

-31.5 

35 

-27.0 

40 

-29.3 

40 

-26.0 

45 

-28.0 

45 

-25.4 

50 

-27.5 

50 

-25.2 

55 

-26.2 

55 

-25.1 

60 

-25.2 

60 

-25.2 

65 

-25.2 

65 

-25.2 

70 

-25.2 

70 

-25.2 

75 

-25.2 

75 

-25.2 

80 
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83 
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87 

-25.2 
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than  25  feet,  cable  effects  may  be  neglected.  However,  the  cable 
capacity  undesirably  influences  the  pulse  formation  and  transmis- 
sion.  For  this  reason,  the  preamplifier  is  located  physically 
close  to  the  BF3  counter.   Besides  its  amplification  effects,  the 
preamplifier  serves  the  purpose  of  impedance  transformation;  that 
is,  to  transform  the  voltage  which  is  developed  across  the  small 
capacity  at  the  input  of  the  preamplifier  into  an  approximately 
equal  voltage  across  the  high  capacity  of  the  output  cable.   The 
preamplifier  used  in  this  work  consisted  of  a  404A  tube  which 
was  d-c  coupled  to  a  417-A  output  cathode  follower.   Capacity 
feedback  could  be  provided  to  the  input  from  either  the  full  out- 
put or  one-third  or  one-tenth  of  the  output.   This  gave  the  pre- 
amplifier gain  a  setting  of  1,  3,  or  10.   Since  most  of  the  time 
a  long  cable  was  used,  the  gain  was  maintained  at  3.   The  entire 
preamplifier  assembly  was  contained  in  the  Detectolab  Model  DA  8 
linear  amplifier.   The  output  of  the  preamplifier  was  coupled 
into  the  differentiating  network  of  the  main  amplifier.   The 
differentiating  network  of  the  amplifier  determines  the  low- 
frequency  cut-off  and  the  clipping  time  of  the  amplifier,  while 
the  integration  network  determines  the  upper  frequency  cut-off 
of  the  amplifier  and  affects  the  rise  and  delay  times  of  the 
pulses. 

The  resolving  time  of  the  proportional  counter,  which  is  the 
minimum  time  that  can  elapse  between  the  start  of  two  pulses  if 
they  are  to  be  counted  as  two  pulses,  is  determined  by  the  rise 
of  the  pulses.   The  resolving  time  of  the  system,  that  is,  the 
counter,  amplifier,  and  scaling  unit,  is  the  result  of  the 
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composite  characteristics  of  all  the  elements  in  it  (19). 

High  voltage  was  fed  directly  to  the  preamplifier  from  the 
high-voltage  power  supply  Model  DVI,  B.  J.  Electronics,  N.  E. 
Dept.  No.  262.   The  preamplifier  furnished  the  signal  to  the 
multi conductor  preamplifier  cable. 

The  gain  of  the  linear  amplifier,  Model  DA8,  B.  J.  Elec- 
tronics, N.  E.  Dept.  No.  376,  consisted  of  a  coarse  attenuator 
with  six  steps,  each  increasing  the  gain  by  a  factor  of  2,  while 
the  fine  attenuator  had  a  continuous  control  and  covered  a  factor 
of  one.   In  order  to  calibrate  this  gain  dial,  a  known  input  sig- 
nal was  introduced  from  the  Radiation  Instruments  Company  pulse 
generator,  Model  1401,  N.  E.  Dept.  No.  1046.   The  output  pulses 

at  each  setting  of  the  gain  dial  were  measured  on  the  screen  of 

Eout 
an  oscilloscope,  Type  Tex.  545.   The  gain,  ,  was  calculated 

Ein 
at  each  setting  of  the  gain  dial  and  the  values  obtained  were 

plotted  against  their  corresponding  dial  setting  In  Fig.  13. 

COUNTER  CHARACTERISTICS 

Proportional  Counter  Action 

Theory.   Consider  an  electron  which  has  been  formed  at  a 
certain  point  within  the  volume  of  the  central  wire  by  an  initial 
ionizing  event.   This  electron  will  drift  towards  the  central 
wire  under  the  influence  of  the  field  due  to  the  applied  voltage. 
As  the  electron  drifts  in,  it  will  make  collisions  with  the  atoms 
and  molecules  of  the  gas  with  which  the  counter  is  filled.   If 
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the  counter  has  cylindrical  geometry,  the  field  varies  inversely 
as  the  first  power  of  the  radius  (19), 

E  =  V0/r  logtrg/rx)  (7) 

where  E  is  the  electric  field  strength,  V0  is  the  applied  volt- 
age, r^  and  ?2   are  ^ne  radii  of  the  wire  and  cathode,  respec- 
tively.  From  equation  (7)  it  can  be  seen  that  the  field  increases 
towards  the  central  wire,  therefore  the  electron  will  acquire 
more  energy  per  unit  path  near  the  central  wire.   Thus  when  the 
electron  gains  enough  energy  in  one  mean  free  path  between  col- 
lisions to  produce  ionization,  the  next  time  it  hits  an  atom  it 
will  ionize  the  atom.   It  is  thus  the  beginning  of  the  process  of 
cumulative  ionization,  which  is  the  basis  of  the  counter  action. 
If  this  collision  resulting  in  ionization  takes  place  at  some 
distance  from  the  central  wire,  the  Initial  electron  plus  the 
new  electron  will  start  on  their  next  free  path  towards  the  cen- 
tral wire.   Since  the  field  increases  nearer  the  wire,  any  elec- 
tron which  has  acquired  enough  energy  to  produce  ionization  by 
collision  at  one  point  will  acquire  enough  energy  to  ionize  in 
subsequent  mean  free  paths. 

For  a  cylindrical  proportional  counter  with  two  concentric 
electrodes,  the  size  of  the  pulse  V  In  volts  appearing  on  the 
central  wire  system  is  given  by  (20) 

V  =  A  ne/C  (8) 

where  A  is  the  gas  multiplication  and  n  is  the  number  of  elec- 
trons in  the  initial  Ionizing  event.   If  V  is  expressed  in  volts, 
then  e  is  the  charge  on  the  electron  in  microcoulombs  and  C  is 
the  distributive  capacity  of  the  central  wire  system  in  micro- 
farads.  The  gas  multiplication  A  is  defined  as  the  number  of 
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additional  electrons  produced  by  each  initial  electron  and  by- 
its  progeny  as  it  travels  from  the  place  where  it  is  originally 
formed  to  the  central  wire.   Thus  A  is  a  measure  of  the  size  of 
the  "avalanche"  which  each  initial  electron  starts,  and  the  aval- 
anches are  assumed  to  be  independent  of  each  other. 

As  the  voltage  V0,  equation  (7),  across  the  counter  is 
slowly  raised,  a  value  is  obtained  at  which  the  multiplication 
starts.   The  lowest  voltage  at  which  this  occurs  is  the  threshold 
voltage  for  a  proportional  counter.   The  physical  significance 
of  this  voltage  lies  in  the  fact  that  it  is  the  minimum  voltage 
at  which  an  electron  traveling  towards  the  central  wire  acquires 
enough  energy  to  ionize  in  its  last  free  path  before  it  reaches 
the  wire.   If  rQ  is  defined  as  the  radius  at  which  the  avalanche 
starts,  r0  will  be  equal  to  r,  at  the  threshold  voltage.   However, 
r0  in  the  region  of  proportional  counting  never  becomes  much 
greater  than  r^  (see  discussion  in  Appendix),  and  therefore  as- 
sumes an  approximately  linear  dependence  on  the  potential.   The 
radius  r0,  at  which  the  avalanche  starts  when  the  voltage  applied 
is  VQ,  is  therefore  related  to  the  threshold  voltage  Vp,  as 
follows  (20): 

r0  =  *l   VVp  <®> 

Thus  rQ  is  defined  in  terms  of  the  directly  measurable  voltages 

and  wire  radius. 

The  multiplication  factor  A,  as  defined  in  equation  (8),  Is 
determined  by  the  number  of  ionizing  collisions  which  the  elec- 
tron makes  as  it  travels  to  the  wire.   Rose  and  Korff  (20)  give 
the  following  general  expression  for  the  gas  multiplication  A  at 
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an  applied  voltage  V0, 

A  =  exp  2  (a  Nm  C  p^(0)  V^|^V0/Vp)*  -  1     (10) 
where  0(0),  the  number  of  slowest  electrons,  Is  to  be  equal  to 
unity  (see  discussion  in  Appendix),  Nm  is  the  number  of  gas  mole- 
cules per  unit  volume,  C  the  capacity  per  unit  length  of  the 
counter,  r^  the  wire  radius,  V  the  threshold  voltage,  and  "a" 
is  the  rate  of  change  of  the  ionization  cross  section  for  ener- 
gies just  above  the  ionization  potential.   Equation  (10)  can  be 
written  as 

A  =  exp  k  V0*  [(Vo/Vp)*  -  ll  (11) 

where  k  =  2  (a  Nm  C  r1  0(O)J*.   This  is  an  adequate  form  for  cal- 
culation of  the  theoretical  values  of  gas  multiplication. 

Experimental  Values  of  Gas  Multiplication.   The  experimental 
values  of  multiplication  for  counter  (B),  Fig.  4,  were  deter- 
mined using  the  circuit  shown  in  Fig.  14. 

The  procedure  for  measuring  the  gas  multiplication  was  as 
follows: 

1.  The  linear  amplifier,  Model  DA  8,  B.  J.  Electronics, 
N.  E.  Dept.  No.  376,  gain  was  reduced  to  a  minimum  and  the  high 
voltage  increased  until  pulses  appeared  on  the  oscilloscope 
screen. 

2.  The  pulse  height-selector  setting  dial  was  increased 
slowly  until  no  pulses  were  visible  on  the  oscilloscope  screen. 
The  PHS  reading  together  with  the  corresponding  high  voltage  and 
the  linear  amplifier  gain  were  noted. 

3.  Voltage  was  reduced  in  steps  of  50  volts,  gain  was 
adjusted  if  necessary,  and  at  each  voltage  setting,  step  2  was 
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repeated. 

4.  The  maximum  pulse  height  corresponding  to  an  arbi- 
trarily selected  point  in  the  plateau  region  of  the  counter  was 
noted. 

5.  The  PHS  readings  obtained  at  the  different  voltage 
settings  were  renormalized  to  the  same  gain.   The  data  so  ob- 
tained was  renormalized  to  a  multiple  of  the  pulse  height  ob- 
tained in  step  4.   Values  of  the  pulse  sizes  thus  obtained  are 
included  in  Table  3  and  plotted  against  the  corresponding  volt- 
ages In  Pig.  15. 

Prom  Fig.  15  it  is  seen  that  the  pulse  size  is  approxi- 
mately constant  between  550  and  750  volts.   Considering  that  the 
counter  was  acting  as  an  ionization  chamber  in  this  region,  that 
is,  the  multiplication  in  this  region  was  unity,  the  absolute 
multiplication  values  were  obtained  by  normalizing  the  arbitrary 
pulse  sizes  obtained  at  different  voltages  by  the  constant  pulse 
size  obtained  in  this  region.   Multiplication  values  thus  ob- 
tained are  listed  in  Table  4  and  plotted  versus  voltage  in 
Pig.  16. 

The  threshold  voltage  Vp  is  the  voltage  corresponding  to 
the  point  where  the  straight-line  portion  of  the  curve,  Pig.  16, 
if  extrapolated  downwards,  would  intersect  the  abscissa.  The 
observed  value  for  Vp  is  855  volts.  The  theoretical  value  of 
Vp  was  calculated  from  equation  (11)  by  considering  the  first 
two  points  lying  on  the  straight- line  portion  of  the  curve, 
Fig.  16,  as  follows. 

If  A^  and  A2  are  the  multiplication  factors  at  these  two 
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BP3  gas  multiplication  data. 

High  voltage 

Normalized  PHS     :     Multiplication 
readings          : 

1750 

1700 

1650 

1600 

1550 

1500 

1450 

1400 

1350 

1300 

1250 

1200 

1150 

1100 

1050 

1000 

950 

900 

850 

800 

750 

700 

650 

600 

550 

0.341                 191.001 
0.267                   149.520 
0.216                  120.960 
0.190                  106.512 
0.148                    63.272 
0.135                   76.600 
0.114                   63.952 
0.083                   46-480 
0.064                    36.176 
0.045                   25.480 
0.031                   17.472 
0.020                   11.200 
0.012                    7.000 
0.009                    5.408 
0.008                     4.536 
0.005                    3.096 
0.005                     2.800 
0.004                     2.598 
0.003                    2.200 
0.003                    2.100 
0.003                    1.971 
0.003                    1.843 
0.002                    1.277 
0.002                    1.199 
0.001                    1.000 

^oo» — 
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points,  V^  and  V2  are  the  corresponding  voltage  values, 
Ax  =  exp  k  VX*  [(Vi/Vp)^  -  1 
A2  =  exp  k  V2i  |"(V2/Vp)*  -  l" 
In  A^  In  Aj 


Therefore 


v^fiVi/Vp)*  -  1]   v2t[(v2/vp)£  -  1 


Using  this  equation,  the  calculated  value  of  Vp  was  795.5 
volts.   Prom  equation  (11)  it  can  be  seen  that  k  is  a  constant 
depending  on  the  nature  and  pressure  of  the  gas  and  on  the  wire 
and  cylinder  radii,   k  also  depends  on  the  energy  distribution 
of  the  secondary  electrons  arising  in  the  avalanche  process.   In 
the  approximation  made  in  deriving  equation  (10),  the  relative 
number  of  slowest  electrons  enters  and  is  adjusted  in  the  param- 
eter a.   Since  no  information  as  to  the  availability  of  data  on 
a  for  BF3  was  present,  the  value  of  k  was  calculated  at  one  volt- 
age V0  by  using  the  theoretical  Vp  value  and  known  multiplication 
factor  A.   Using  this  method,  the  value  of  k,  calculated  to  be 
0.034089,  was  used  to  calculate  the  theoretical  values  of  multi- 
plication at  different  voltages.   Theoretical  values  of  multipli- 
cation obtained  as  such  are  listed  in  Table  5  and  plotted  in 
Fig.  16. 

It  is  to  be  noted  that  equation  (11)  predicts  that  values 
of  log  A  corresponding  to  different  V0  values,  contrary  to  Fig. 
16,  should  lie  on  a  straight  line  of  slope  (d  log  A/dV0). 

This  disagreement  is  apparent,  however,  only  when  one  takes 
into  consideration  the  assumption  on  which  the  derivation  of 
equation  (10)  Is  based.   The  assumption  of  interest  in  the  present 
connection  involves  neglecting  of  fluctuations  in  energy  loss,  as 
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It  is  assumed  that  at  a  given  cylinder  voltage  there  exists  a 
point  at  a  distance  r0  from  the  wire  axis  at  which  the  avalanches 
start,  and  that  this  distance  is  the  same  for  all  electrons  ini- 
tiating such  avalanches.   The  value  of  r0  which  was  used  in  ob- 
taining equation  (10)  was,  in  fact,  the  average  starting  point 
for  sustained  ionization.   While  this  assumption  would  be  ex- 
pected to  be  valid  for  large  voltages,  for  which  the  mean  r0  is 
considerably  greater  than  the  wire  radius  r1#  it  should  also  be 
expected  that  it  will  break  down  at  low  voltages.   This  is  im- 
mediately evident  when  one  considers  the  effect  of  a  distribution 
in  starting  points  r0  which  would  naturally  enter  as  en  essential 
component  of  the  avalanche  model  when  fluctuations  are  taken 
into  account. 

At  low  voltages  the  gas  multiplication  is,  at  first,  caused 
by  these  few  electrons  which  have  large  rQ.   As  the  voltage  is 
increased,  the  r0  distribution  function  shifts  outward  from  the 
wire  until  it  Is  virtually  outside  the  point  v-^   so  that  practi- 
cally all  the  electrons  contribute  to  the  multiplication.   There 
fore  the  amount  of  charge  collected  when  gas  multiplication  sets 
in,  and  hence  the  multiplication  itself,  would  increase  more 
slowly  than  If  all  avalanches  started  at  the  same  point  and  a 
curve  of  the  sort  obtained  (d  log  A/dV  or  dA/dV  ■  0  for  A  =  1) 
would  be  expected. 

Hence  It  can  be  stated  that  the  form  of  the  multiplication 
curve  near  the  threshold  is  to  a  certain  extent  a  picture  of  the 
r0  distribution  function  for  r0  larger  than  the  average  value. 
Therefore  at  higher  voltages  when  all  but  a  small  portion  of  r0 
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distribution  function  lies  outside  r^,  almost  all  the  primary 
electrons  are  effective  in  initiating  avalanches.   Thus  the  as- 
sumption of  a  fixed  r0,  while  not  rigorously  fulfilled,  gives  an 
asymptotic  value  of  A  in  equation  (11)  for  very  large  voltages. 

Sensitivity  of  Boron  Trifluoride  Counters 

Theory.   The  B10  (n,  oc )  reaction  has  been  studied  quite  ex- 
tensively by  Brower,  Bretscher,  and  Yilhent  (8).  The  reaction  is 
exothermic,  with  an  energy  release  of  2.78  mev.   The  ground  state 
of  the  Li  product  nucleus  may  be  formed  directly,  with  the  en- 
tire energy  being  shared  by  the  Li7  and  the  alpha  particle. 
Alternatively,  an  intermediate  excited  state  of  Li  may  be  formed, 
followed  by  the  emission  of  a  0.48  mev  gamma.   The  probability 
of  reaching  the  ground  state  directly  is  only  0.07  for  reactions 
induced  by  thermal  neutrons.   It  rises  to  a  maximum  of  about  0.7 
for  1.8  mev  neutrons,  followed  by  a  drop  to  less  than  0.5  for 
2.5  mev  neutrons  (18). 

For  the  capture  of  a  neutron  of  energy  E,  the  kinetic  energy, 
Eke,  of  either  E  +  2.30  mev  or  E  +  2.78  mev  will  be  shared  by  the 
alpha  particle  and  the  lithium  recoil  nucleus.   If  the  captured 
neutron  is  thermalized,  E  is  approximately  zero,  and  the  energies 
of  the  lithium  nucleus  and  the  alpha  particle  in  this  case  are 
(19): 

Eke   Mo^         (2.30)(4.00) 

ELi  ■  -  a  0.88  mev 

Moc  +  MLi        4.00  +  7.02 

and 
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Eke  MLi 
Eo^=  =  1.47  mev 

MoC  +  MLi 
where  McXL  and  M^  are  the  mass  of  c<  -particle  and  Li-nucleus, 

respectively. 

The  energy  dependence  of  the  cross  section  q~   for  the  B 

(n,  o(  )  reaction  was  found  to  vary  inversely  with  the  velocity 

to  about  100  ev,  and  therefore  can  be  represented  mathematically 

as 

°o  vo 

cr« (12) 

v 
where  CTQ   is  the  cross  section  of  the  reaction  of  B10  atoms  with 
the  thermal  neutrons  at  room  temperature,  having  a  value  of 
4,010  barns,  and  v0  is  the  velocity  of  thermal  neutrons  at  293° 
K  with  a  value  of  2.2  x  105  cm  per  second. 

The  reaction  rate  in  such  a  detector  can  be  obtained  by  (19) 

n  v      _    <-. 

r  =  _  =  n  vl_  =  0Z_ 

where  n  is  the  number  of  neutrons  per  cubic  centimeter,  R  is  the 
reaction  rate  per  cubic  centimeter,  0   and  Zl  are  the  neutron 
flux  and  the  total  cross  section  of  B10  (n,  <?0  reaction.   The 
contribution  dR  to  the  reaction  rate  by  neutrons  between  k  and 
E  +  dE  in  the  volume  dV  is 

dR  =  N(X,  Y,  Z)  0(E,    X,  Y,  Z)  (T(E)  dEdV 
where  N(X,  Y,  Z)  is  the  number  of  B10  atoms  per  unit  volume  at 
point  X,  Y,  Z,  0~(E)  is  the  cross  section  at  energy  E,  and 
0(E,   X,  Y,  Z)  is  the  neutron  flux  per  unit  energy  interval  at 
the  point  X,  Y,  Z. 
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The  total  reaction  rate  throughout  the  entire  detector  could 
be  represented  by  the  integral  over  the  detector  volume  and  the 
complete  neutron-energy  spectrum,  and  can  be  written  as 

Reactions  per  second  =       f   N(X,  Y,  Z)  (J  (E) 


1  vol  'o 


0{E,    X,  Y,  Z)  dVdE         (13) 
If  the  flux  can  be  considered  independent  of  the  position  in  the 
detector,  and  if  the  total  number  of  B^  nuclei  can  be  repre- 
sented by  Nrp,  equation  (13)  becomes 

Reactions  per  second  =  NT  I  <T~(E)  0{E)   dE        (14) 

For  a  neutron  spectrum  between  0  and  100  ev,  equation  (12)  can 

be  applied  and  equation  (14)  becomes  (19) 

^100  ev 
Reactions  per  second  =  Ny  0~o   v0  )       n(E)  dE 

=  NT  0"0  vQ  n  (15) 

where  n  is  the  number  of  neutrons  per  unit-volume,  irrespective 
of  energy  up  to  100  ev.   If  P(E)  represents  the  probability  that 
a  B1  (n,  o<  )   reaction  caused  by  a  neutron  of  energy  E  will 
register  a  count,  then  from  equation  (15)  the  expression  for  the 
counting  rate  is 

Counts/sec  =  NV  0~0   vQ  j   P(E)  n(E)dE  (16) 

where  NT  ■  NV,  that  is,  the  product  of  the  density  of  B10P3 
atoms  and  the  volume  of  the  detector. 

In  most  cases,  the  counter  is  used  with  a  counting  system 
that  requires  a  certain  minimum  pulse  size  for  it  to  be  triggered. 


51 


This  means  that  the  total  energy  of  the  alpha  particle  and  the 
lithium  nucleus  dissipated  in  the  chamber  should  exceed  a  certain 
threshold  value  Et.   Since  in  practice  the  chamber  dimensions  are 
usually  larger  than  the  particle  ranges  and  the  threshold  energy 
Et  is  much  less  than  the  total  energy  available,  the  probability 
function  P(E)  can  be  taken  as  1  to  a  good  approximation. 

If  the  sensitivity  of  a  counter  is  defined  as  the  counts  per 
second  per  unit  neutron  flux,  the  expression  for  sensitivity  as 
given  by  equation  (16)  is 

Sensitivity  =  NV  0~Q   (Vq/v)  (17) 

For  Maxwell- Bo It zman  distribution,  this  becomes 

4,010  x  10~24  NV 


Sensitivity  ■ 


1.128 


where  v  =  1.128  vQ,  N  is  the  total  number  of  B10  atoms  per  cubic 
centimeter,  and  V  is  the  sensitive  volume  of  the  counter.  Equa- 
tion (17)  is  valid  only  if  the  flux  is  constant  throughout  the 
counter  volume. 

Experimental  Determination  of  Sensitivity.   The  experimental 
value  of  the  sensitivity  was  determined  by  placing  the  counter  at 
different  positions  in  the  Kansas  State  University  graphite  sub- 
critical  assembly.   Positions  of  the  slots  and  source  are  shown 
in  Pig.  17.   This  pile  consists  of  a  rectangular  parallelepiped 
68  inches  square  and  100  inches  high,  which  rests  on  a  concrete 
foundation.   Machined  reactor  grade  graphite  blocks,  about  four 
inches  in  cross  section  and  of  various  lengths,  were  used  to 
construct  the  column.   Some  of  the  graphite  blocks,  as  shown  in 
Pig.  17,  were  drilled  to  a  diameter  of  1.75  inches  along  the 
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Fig.  17.   Kansas  State  University  subcritical  assembly. 
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length  of  the  block  to  accommodate  fuel  slots.   Each  of  these 
fuel  slots  contained  graphite  cylinders  measuring  1.625  inches 
in  diameter  and  22.68  inches  in  length.   These  cylinders,  called 
plugs,  are  made  of  the  same  material  as  the  graphite  blocks. 

For  the  standardization  of  counter  (B),  Pu-Be  neutron  source 
#3761,  containing  about  16  gms  of  Pu  and  emitting  approximately 
1.71  x  10°  neutrons  per  second,  was  placed  in  the  center  of  fuel 
port  1,  Pig.  17.   Accurate  determination  of  the  flux  distribu- 
tion in  the  graphite  pile  using  indium  foils  and  small  commercial 
BF3  counters  was  carried  out  by  Foulke  (12),  using  a  Pu-Be  source 
emitting  1.58  x  10"  neutrons  per  second.  A  conversion  factor  of 
1.0823  was  therefore  used  to  correct  for  the  difference  in  the 
source  strengths. 

The  counter  was  connected  directly  to  a  Nuclear  Chicago  pre- 
amplifier, Model  1062,  N.  E.  Dept.  No.  96,  output  of  which  was 
fed  to  a  Baird  atomic  scaling  unit,  Model  132,  N.  E.  Dept.  No. 
145.   A  regulated  high-voltage  supply,  John  Fluke  Model  400BDA, 
N.  E.  Dept.  No.  179,  was  used  as  the  high-voltage  source  for 
the  counter. 

The  measurements  were  run  so  as  to  give  10,000  counts  in 
each  position.   The  count  rates  were  corrected  for  background 
and  listed  in  column  4  of  Table  6.   The  experimental  values  of 
sensitivity  at  each  position  of  the  pile  were  independently  cal- 
culated.  Values  obtained  are  given  in  column  5  of  Table  6,  and 
plotted  against  distance  from  source  in  Fig.  18.   From  this  figure 
it  can  be  seen  that  the  sensitivity  increases  with  the  increase 
in  distance  up  to  about  1.5  diffusion  lengths  and  then  remains 
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constant.  Experiments  done  In  this  department  on  the  graphite 
assembly  show  that  the  Cadmium  Ratio,  defined  as  the  ratio  of 
total  flux  to  epicadmium  flux,  rises  sharply  at  approximately 
1.5  diffusion  lengths  from  the  source  which,  for  graphite,  cor- 
responds to  about  75  cm.   Thus  the  increase  of  sensitivity  with 
the  increase  in  distance  can  be  attributed  to  the  fact  that  at 
smaller  distances  neutrons  are  not  thermalized  enough  and  that 
the  cross  section  of  B10  (n,  oC  )  reaction  is  correspondingly 
lower . 

Equation  (17)  was  used  to  calculate  the  sensitivity  of  the 
counter  as  such: 

Pressure  of  BP3  gas  inside  the  counter  =  29.9  cm  Hg 

Sensitive  length  =  5  inches 

Diameter  of  the  cylinder  =0.94  inch 

Sensitive  volume  =  irR2h  =  56.8079  cm3 
The  sensitive  volume  has  been  calculated  here  on  the  assump- 
tion that  it  is  a  rectangular  cylinder.   Since  the  central  wire 
was  supported  at  one  end  by  a  thicker  metal  rod,  the  field 
strength  at  the  surface  of  the  wire  is  weaker  near  this  end  and 
reaches  its  normal  value  only  at  some  distance  from  this  end 
(19).   The  field  near  the  end  is  therefore  different  from  the 
field  near  the  center  and  the  lines  of  forces  are  not  radial  so 
that  the  sensitive  volume  is  not  a  rectangular  cylinder  but  has  a 
conical  shape  near  this  end.   Since  the  boron  trifluoride  gas  was 
96  per  cent  enriched  in  B10  component,  the  density  N  of  the  B10 
atoms  in  the  counter  was: 
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(6.02  x  1023)  (atoms/gin  atom)  (29.96)  (pressure,  cm  Hg)(0.96) 

(atoms  B^O/total  atoms) 

N  = 

(76) (atmospheric  pressure,  cm  Kg)(22.4  x  103) (cm3/gm  mole) 
=  9.8954  x  1018  atoms  per  cm3 
Substituting  the  value  of  V  and  N  in  equation  (17), 

(4010  x  10~24  NV) 


Sensitivity  = 


1.128 
4010  x  10~24(0.9895  x  1019) ( 56.8079) 


1.128 
=  1.998  counts  per  unit  neutron  flux  per  second 
Comparison  of  this  result  with  the  experimental  values  in 
Table  6  shows  that  the  latter  are  between  4  and  5  per  cent  of  the 
theoretical  values.   This  discrepancy  can  be  attributed  to  sev- 
eral factors,  among  which  are  the  inaccuracy  in  determining  the 
sensitive  volume  (19),  the  thickness  of  the  anode  wire  (19), 
the  gas  impurities  (4),  and  the  electronic  circuitry. 

Plateau  and  Integral  Bias  Curves 

Some  typical  curves  of  counting  rate  versus  voltage  applied 
to  the  counter  are  shown  in  Figs.  19  and  20.   The  neutron  source 
used  to  obtain  these  data  were  a  Pu-Be  source  moderated  In  par- 
affin.  A  linear  amplifier,  type  DA8,  B.  J.  Electronics,  N.  E. 
Dept.  No.  364,  in  conjunction  with  a  scaler  DS9,  B.  J.  Electron- 
ics, N.  E.  Dept.  No.  144,  formed  the  electronic  circuit.   The 
counter  was  connected  to  the  preamplifier  through  a  coaxial  cable 
type  RG11-U.   Figures  19  and  20  show  the  effect  of  amplifier  gain 
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on  the  plateau  characteristics.   As  shown  in  Pig.  19,  the  plateau 
is  about  250  volts  at  a  gain  of  8;  however,  it  becomes  shorter 
and  steeper  with  higher  gain. 

An  integral  bias  curve  with  1475  volts  applied  to  the  count- 
er is  shown  in  Pig.  21.   The  relatively  broad  distribution  in  the 
count  rate  is  attributed  to  the  partial  ionization  loss  by  the 
disintegration  products  which  strike  the  wall  of  the  counter. 
The  geometry  is  favorable  for  the  production  of  these  small  pulses 
since  the  range  of  the  disintegration  product  is  of  the  order  of 
half  the  counter  radius.   If  the  ratio  of  range  to  the  counter 
radius  were  decreased  by  using  a  larger  diameter  or  higher  fill- 
ing, the  effect  would  be  considerably  reduced.   However,  at  the 
same  time  the  sensitivity  to  the  gamma  radiation  would  be  in- 
creased, an  important  consideration  for  reactor  application. 

Effect  of  Gamma  Rays  on  the  Counter  Performance 

In  order  to  determine  the  relative  efficiency  of  the  counter 
for  gamma  rays  as  compared  to  its  efficiency  for  neutrons,  and, 
further,  to  determine  the  effect  of  the  amplifier  gain  on  this 
relative  efficiency,  the  following  experiment  was  performed. 
1.  The  counter  was  placed  at  a  fixed  distance  from 
Nuclear  Chicago  moisture  gauge,  Model  P-21,  N.  E.  Dept.  No.  136, 
containing  five -mi 111 curie  Ra-Be  source  surrounded  by  paraffin 
to  provide  a  beam  of  thermal  neutrons.   Counts  were  taken  at  sev- 
eral 5-minute  intervals  at  a  certain  setting  on  the  amplifier. 
The  average  count  rate  corresponding  to  that  particular  gain 
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setting  was  determined,  corrected  for  background,  and  recorded 
in  column  3  of  Table  7.   Background  and  background  plus  neutron 
count  rates  are  plotted  versus  gain  dial  setting  in  Fig.  22. 

2.  The  dose  rate  received  by  the  counter  from  the 
source  was  determined  by  placing  a  slow-fast  survey  meter,  Model 
RCL  20806,  N.  E.  Dept.  No.  10178,  at  the  same  distance  from  the 
source.   The  dose  rate  obtained  in  this  case  was  0.0046  mr  per 
hour. 

3.  Procedure  adopted  in  step  1  was  repeated  for  gamma 
rays.   Source  used  in  this  experiment  was  a  Cs-137  one,  and 
placed  at  four  inches  from  the  counter.   The  average  count  rate 
versus  gain  dial  setting  on  the  amplifier  is  plotted  in  Fig.  22, 
and  the  values  corrected  for  background  are  listed  in  column  2 
of  Table  7. 

4.  The  dose  rate  received  by  the  counter  from  the  gamma 
ray  source  was  determined  by  a  GM  survey  meter,  SM-136,  serial 
1408,  N.  E.  Dept.  No.  10156.   The  value  obtained  was  6.5307  mr 
per  hour,  at  four  inches  from  the  source. 

5.  Count  rate  values  for  neutrons  and  gamma  rays  were 
normalized  to  the  same  source  strength  and  listed  in  columns  5 
and  4  of  Table  7. 

6.  The  normalized  gamma  counts  were  divided  by  the  cor- 
responding normalized  neutron  readings  to  determine  the  relative 
counter  efficiencies.   Results  are  listed  in  column  6  of  Table  7 
and  plotted  in  Fig.  23  against  the  corresponding  gain  settings 

in  column  1.   Figure  23  shows  that  the  relative  counter  efficiency 
for  gamma  rays  increases  with  the  increase  in  gain.   This  can  be 
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attributed  to  the  fact  that  the  average  pulse  height  produced 
by  a  gamma  ray  is  of  the  order  of  one-hundredth  of  the  pulse 
height  produced  by  a  thermal  neutron.   At  low  gains  the  scaler 
cannot  see  the  0    pulses,  but  at  higher  gain  the  gamma  pulses 
are  amplified  enough  to  trigger  the  scaler,  and  therefore  a 
higher  number  of  counts  due  to  gammas  is  recorded. 


USE  OP  THE  COUNTER  IN  DETERMINING  THE  DIFFUSION 
LENGTH  OF  THERMAL  NEUTRONS  IN  THE  KANSAS  STATE 
UNIVERSITY  GRAPHITE  MODERATED 
SUBCRITICAL  REACTOR 


Theory 

The  relationship  between  the  thermal  diffusion  length  and 
the  thermal  flux  distribution  in  a  homogeneous  moderating  medium 
according  to  a  raonoenergetic  model  is  given  by  (11), 

V     0   -  -s  =  0  (18) 

L2 

where  0   is  the  thermal  neutron  flux,  and  L  is  the  diffusion 

length.  Equation  (18)  is  valid  for  the  thermal  flux  in  regions 

far  enough  from  the  source  so  that  there  is  no  production  of 

thermal  neutrons,  and  a  few  mean  free  paths  from  the  boundaries 

of  the  pile.   The  solution  of  equation  (18)  is  given  by  (11), 

OO  <yO 

-v—  ^--         mrrx     ntry 

0  -   2L_  Z-  A^i  cos  cos  sinh  }f   (c  -  z)   (19) 

m=l  n=l         a       b 

where  a  and  b  are  the  effective  x  and  y  dimensions  of  the  pile, 

respectively,  c  is  the  effective  pile  height  measured  from  the 
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source  plane,  z  is  the  distance  from  the  source  plane,  l/J^  is 
the  relaxation  length  of  thermal  neutrons  of  the  mn  th  harmonic 
in  the  z-direction  given  by, 

^.o    1     mtr  o    nrr  o 

C  "  75  +  «-»    +  <->2  (2o» 

Ifi  a       b 

and  Aj^  is  a  function  of  the  assumed  neutron  source  boundary 
condition  in  the  pile.   Considering  a  point  thermal  neutron 
source  emitting  S  neutrons  per  second,  and  represented  by  a  Dirac 
delta  function,  the  center  of  source  plane  being  taken  as  the 
center  of  the  coordinate  system,  A^,  is  given  by  (11) 

2   S 

Aran  =  ~7Ty T7 (21) 

abD  Jmn  cosh^mn  c 
where  D  is   the  diffusion  coefficient  of  thermal  neutrons  in  the 

lattice. 

Combining  equations    (19)    and   (21)    and  rearranging  terms, 

the   solution  of  equation   (18)    can  be  written   in  the  form  (11), 


2S     ^  2?  1  mrrx 

0{x,  y,    z)  =  /__  C.  "T —  cos  

abD     m=l,  3,5,...      n=l,3,5, ...     "mn  a 


linn^™,,    (c  -   z) 


nrry     lU.  v  xan 

cos  (22) 

b       cosh  ^     • 

By  rewriting  sinh  /^  (c  -  z)  in  exponential  form  and  rearrang- 
ing terms,  equation  (22)  can  be  written  as  (15), 

2  S      ^  eS  1  mirx 

fix,  y,   z)  =  2-  2- 


cos 


abD     m=l,  3, 5,...      n=l,3,5,...     ^nm 

mn^c-z) 


niry    -  y  z 

cos  e       n*11     x 

b 


1  +  e"2  ^nm  c 


(23) 


By  expanding  equation   (23),    the  expression  for  0(x,    y,    z)    can  be 
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written  as, 


0(x,  y,  z)  = 


2  S 


11 


abD 
7 

*Z    Z. 

m=l  n=l 
(except  1,1) 


1       rrx 

— —  cos  COS 

0-n     a      b 


—  e-^ll25  (FMECF) 


mtrx     ntry  _y 


y- 

,/  mn 


cos  cos  e 

a      b 


mn' 


.  -  e"2^inntc-z) 


TV** 


mn 


(24) 


where  FMECF  is  the  Fundamental  Mode  End  Correction  Factor  and 

•m  (15), 
1  -  e-^111 


is  written  in  the  form  (15), 
FMECF  = 


1  +   e-2JH  c 
Equation   (24)    can  be  rearranged  as   follows: 


(25) 


0(x,    y,    z) 


2  S 
abD 

x     1  + 


r- 


cos 


*x  try        V  > 

—  cos  —  e      11     (FMECF) 


11 


/ne^lZ   (1+  e"2^lC) 


cos 


—  cos  —  1  -   e        **  / 

a  b    ^ 


^jL.       JL  mtrx  nrnry 

/         f"  cos  cos  

m=l     n=l  a  b 

(except   1,1) 

e-^mnz    A    „   e-2^mn(c-z) 


mn 


1  +  e-^mn  c 


(26) 


The   above   equation  can  be  put   in  the  form   (15), 


2  S 

0(x,y,z)   ■  

abD 


t; 


rrx  try        Y 

cos  —  cos  —  e"    HZ    (FMECF)  (CEHCFU 


(27) 


11     a      b  j 

where  CEHCF  is  the  Composite  End  and  Harmonic  Correction  Factor 
and  is  given  by  (15), 
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CEHCP  = 


1    + 


TTX 
COS    —   COS 

a  b 


!!(!  .   e-2^ii(c-.Z)j 


7  7  e"^nm  z     /i   .   Q-2^Wc-z) 


:^% 


m=l     n=l  m         \    1  +  e"'  'a*  c 

(except   1,1) 

mirx           ntry 
cos  cos  

a      b 


(28) 


The  thermal  diffusion  length  is  related  to  the  relaxation  length 
by  (11), 

The  diffusion  length  L  is  calculated  from  the  flux  data 
taken  along  the  z-axis  of  the  pile  by  using  equations  (24)  to 
(28).   The  general  procedure  for  calculating  L  is  as  follows. 

1.  A  provisional  value  of  ]f..   is  obtained  from  the 
linear  slope  of  the  In  0(z)    against  z  curve  for  points  not  too 
near  the  neutron  source  or  the  top  of  the  subcritical. 

2.  With  this  value  of  y^,  the  If  ma,    FMECF,  and  CEHCF 
are  calculated,  using  equations  (20),  (25),  and  (28),  respec- 
tively. 

3.  A  new  value  of  ^.l  is  then  computed  from  a  least 
squares  fit  of  flux  data  corrected  by  equation  (27). 

4.  With  this  new  value  of  0  \\>    the  process  is  repeated 
until  the  value  of  If y±   calculated  from  equation  (27)  coincides 
with  the  value  used  to  calculate  CEHCP  and  FMECF. 

5.  The  value  of  If  yl   thus  obtained  is  used  to  calculate 
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the  diffusion  length  using  equation  (29). 

To  eliminate  tedious  computations,  the  above  procedure  has 
been  programmed  through  the  7th  harmonics  for  the  IBM  650  digi- 
tal computer  by  Kaiser  (15).   This  program  was  used  to  calculate 
the  diffusion  length  using  data  obtained  by  the  author. 

Experimental 

Experimental  determination  of  the  diffusion  length  L  was 
carried  out  as  follows. 

Five  one-curie  (Pu-Be)  sources  encapsulated  in  tantalum  and 
stainless  steel,  emitting  approximately  8.0  x  106  neutrons  per 
second,  were  placed  in  the  Kansas  State  University  subcritical 
assembly. 

The  BP3  counter  was  connected  directly  to  a  Nuclear  Chicago 
preamplifier,  Model  1062,  N.  E.  Dept.  No.  96,  and  was  placed  in 
the  center  of  each  fuel  port,  as  shown  in  Fig.  17.   Count  rates 
C(z)  at  the  different  z  positions  were  found  by  determining  the 
counting  time  for  10,000  preset  counts,  and  were  recorded  in 
column  2  of  Table  8. 

In  C(z)  values  were  plotted  against  z  in  Fig.  24  and  the 
value  of  0^   obtained  was  0.06810  cm"1.   This  value  of  /-q  was 
used  in  the  last  section. 

The  extrapolated  dimensions  of  the  pile  a  and  b  were  deter- 
mined by 

a  =  a0  +  2(0.7l)At 
where  aQ  is  the  actual  dimension  of  the  pile  and  A t  is  the 
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Table  8.  Data  for  diffusion  length  measurement. 


Distance 
from 
source, 
z  inch 

:   Counts  per 
:  minute 

: 

:   CEHCF    : 

•  • 

•  • 

:          : 

•  • 

•  • 

FMECF    : 

I  Corrected 
1  counts  per 
:  minute 

20 

18976.40 

•  •  • 

•  •  • 

... 

28 

11687.40 

•  •  • 

•  •  • 

... 

36 

6586.40 

1.086107 

0.9998805 

6064.9507 

44 

3534.40 

1.046538 

0.9995840 

3378.6340 

52 

1856.40 

1.025585 

0.9985539 

1812.7089 

60 

986.62 

1.014250 

0.9949738 

977.6714 

68 

587.66 

1.008077 

0.9825354 

593.3133 

76 

259.65 

1.004756 

0.9393219 

275.1142 

84 

110.65 

1.003091 

0.7891815 

140 . 1554 

92 

18.71 

•  •  • 

•  •  • 

•  •  • 

L  =  57.322  cm  +  ?   (One  set  of  such  data  only  was  taken.) 


transport  mean  free  path.  With  a0  =  b0  ■  172.448  cm  and 
7^  =  2.46,  the  extrapolated  dimensions  are  175.941  cm. 

Corrected  values  of  the  count  rate  are  listed  in  column  5 
of  Table  8.   Values  of  the  CEHCF  and  FMECH  are  listed  in  columns 
3  and  4  of  the  same  table.   Count  rates  obtained  at  points  very 
near  and  very  far  from  the  sources  were  neglected  in  the  IBM  650 
calculations.   Points  lying  on  the  straight-line  portion  of 
Fig.  24  were  the  only  ones  considered.   The  L  value  obtained  was 
57.32  cm  based  on  one  single  set  of  data. 

Kaiser  (15)  reported  an  L  value  of  52.72  cm,  based  on  a  more 
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thorough  experimental  investigation  on  the  same  graphite  assembly 
using  small  BP3  Nuclear  Chicago  counters  and  indium  foils. 

An  L  value  much  closer  to  the  figure  reported  by  Kaiser  (15) 
could  have  been  obtained  had  more  data  been  taken.   Since  the 
purpose  of  the  experiment  was  only  to  demonstrate  the  quality  of 
the  counter,  no  further  attempts  were  made. 

CONCLUSIONS 

The  following  conclusions  are  drawn  from  the  results  of  the 
present  investigation. 

1.  It  is  possible  to  build  up  boron  trifluoride  count- 
ers with  good  characteristics  using  a  relatively  inexpensive 
vacuum  system. 

2.  Counters  should  be  cleaned  very  thoroughly  before 
being  sealed  onto  the  vacuum  system  and  should  then  be  heated  to 
about  250°  C  for  12  hours,  then  to  400°  C  for  two  hours  while 
maintaining  the  pressure  below  10"4  mm  Hg  to  drive  off  occluded 
gases. 

3.  Glass  stopcocks  and  stopcock  grease  in  particular 
should  be  completely  avoided  in  vacuum  systems  used  for  filling 
BP3  counters. 

4.  Experimental  sensitivity  of  the  counter  prepared 
was  much  less  than  its  theoretically  calculated  sensitivity, 
probably  due  to  the  eccentricities  in  anode  wires,  incomplete 
cleaning  of  counter,  and/or  inaccurate  evaluation  of  the  counter 
sensitive  volume. 
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5.  Efficiency  of  the  prepared  counter  for  gamma  rays 
was  of  the  order  of  10"*4  of  its  efficiency  for  neutrons.  In- 
crease in  amplifier  gain  raised  the  counter  efficiency  for  gammas 
but  was  always  below  5  x  10   of  the  efficiency  for  neutrons. 

6.  Determination  of  diffusion  lengths  of  thermal  neu- 
trons in  graphite  using  the  prepared  counter  yielded,  based  on 
the  one  set  of  data  obtained,  a  value  reasonably  close  to  that 
determined  by  indium-foils  irradiation  technique.   Better  values 
could  have  been  achieved  by  improving  counting  statistics  and 
by  taking  the  depression  in  neutron  flux  caused  by  the  presence 
of  this  relatively  large  counter  in  consideration. 
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APPENDIX1 

Determination  of  Gas  Multiplication  Factor 

The  expression  for  the  multiplication  factor  in  terms  of  the 
counter  parameters  has  been  deduced  by  Rose  and  Korff  (20).   The 
equation  for  the  gas  multiplication  was  obtained  under  the  fol- 
lowing assumptions. 

1.  Secondary  electron  emission  by  positive  ions  at  the 
cathode  is  negligible. 

2.  Recombination  and  electron  attachment  to  neutral 
molecules  are  negligible.   The  recombination  process  is  unimpor- 
tant under  high  fields  and  in  low  pressure  counters.   Electron 
attachment  for  some  gases  is  quite  appreciable  (16). 

3.  Photoelectric  effect  at  the  cathode  is  negligible. 
For  polyatomic  molecules  this  assumption  holds  fairly  well. 

4.  The  fluctuation  In  avalanche  process  can  be  neg- 
lected.  Considering  that  a  single  electron  arriving  from  the 
initial  ionization  produces  further  electrons,  and  denoting  the 
number  of  such  electrons  at  a  distance  r  from  the  center  wire  by 
N(r),  the  Integral  equation  for  the  Increase  of  N(r)  when  r  de- 
creases by  dr  can  be  written  as  (20): 

femax  n(<£,  r) 
-dN  ■  N(p)  dr  )      d6  (i) 


'0 


K£) 


^oted  from  Rose,  M.  E.  and  S.  A.  Korff.   Phys.  Rev.  59: 
850.   1941. 
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Here  6"  denotes  the  excess  of  the  electron  energy  over  the  (aver- 
age) ionization  potential  of  the  gas  in  the  counter.   For  the 
sake  of  brevity,  the  excess  energy  £  is  referred  to  as  "the 
energy".  In  equation  (1)  n(6  ,  r)d<£  ia  the  fraction  of  elec- 
trons at  r  with  energy  between  6  and  C   +  d£,  that  is,  the 
energy  distribution  normalized  to  unity.   1( £  )  is  the  mean  free 
path  for  the  ionization  and  <=■  max  is  the  maximum  energy  of  the 
electron  at  r.   The  integral  in  (1)  is  thus  simply  the  average 
number  of  ionizing  collisions  per  unit  distance  and  is  denoted 
by  K(r),  so  that  from  equation  (1) 

N(r)  =  exp      K(r)  dr  (2) 

'r 

In  equation  (2),  r0  is  introduced  as  the  distance  from  the  center 
of  the  wire  at  which  the  avalanche-producing-ionization  starts. 

The  amplification  factor  A  is  then  given  by  N(ri)  where 
r-^  is  the  radius  of  the  wire, 

A   =  N(rx)  =  exp  J    K(r)  dr  (3) 

It  is  seen,  as  may  have  been  anticipated,  that  the  main  task 
of  the  theory  is  the  calculation  of  the  average  number  of  ioniz- 
ing collisions  per  unit  distance,  K(r),  and  the  spatial  extent  of 
the  ionization  region  r0.   Considering  the  former  quantity,  it  is 
clear  that  in  principle  a  knowledge  of  the  energy  distribution  at 
each  point  r  is  involved,  although  all  that  is  eventually  re- 
quired is  the  average  value  of  a  known  function  of  the  energy, 
viz: 

K(r)  =  [l/i(£)]AV  (4) 
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In  order  to  specify  the  energy  dependence  of  the  mean  free 
path,  it  should  be  recognized  that  the  electrons  attain  an  energy 
of  the  order  of  the  ionization  potential  (^15  volt}  at  a  point 
very  near  the  wire  and  that  the  energy  gain  in  the  region  of 
ionization  is  sufficiently  slow  so  that  the  electron  energy  at 
the  wire  is  only  a  few  times  the  ionization  potential.   The 
energy  region  pertinent  to  our  consideration  is  such  that  the 
ionization  cross  section  increases  linearly  with  energy.   Hence 
writing  for  the  cross  section 

CT(E)  =  a£  (5) 

where  a  is  a  constant  depending  on  the  nature  of  the  gas,  and 
for  the  mean  free  path 

i(€)  =  l/aKm£  (5a) 

where  Nm  is  the  number  of  molecules  per  unit  volume,  equation 
(4)  takes  the  form, 

K(r)  =  aNm£AV  (6) 

In  order  to  obtain  the  average  energy  ^Ay»  the  manner  in  which 
the  energy  distribution  arises  should  be  considered.   The  dis- 
tribution, according  to  the  present  scheme  (no  fluctuations), 
arises  from  the  fact  that  the  electrons  at  r  have  different 
points  of  origin  r'  and  the  energy  at  r  is  given  by  the  average 
energy  gain  in  traversing  the  distances  from  r'  to  r.   The 
average  energy  gain  is  determined  by 

d£ 

=  (C  Vo/r)  -  /(6  )  (7) 

dr 

Equation  (7)  determines  £  for  those  electrons  at  r  which  have 

started  from  r',  as  a  function  of  r  and  r'  subject  to  the 
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initial  condition 

<£(r  =  r')  =  0 
Here  %(£■)    is  the  energy  loss  per  unit  distance  due  to  all  en- 
ergetically possible  inelastic  processes,  V0  is  the  voltage 
across  the  counter,  and  l/C  =  logt^/r^)  so  that  $C  is  the 
capacity  per  unit  length  of  the  counter  with  r2  the  cylinder 
(cathode)  radius. 

According  to  the  above,  the  number  of  electrons  at  r  having 
an  energy  between  (£.   and  €.   +  d  £     is  just  the  number  which  was 
formed  at  a  distance  between  r'  and  r'  +  dr',  where 
dr!  =  (dr'/dOcl£  and  r'  is  given  in  terms  of  £  and  r  by  equa- 
tion (7).   Thus  the  energy  distribution  at  r  depends  on  the  dis- 
tribution at  r',  and  is  accordingly  given  by  the  following  in- 
tegral equation: 

re  max  n  (  £  f ,  r '  ) 
N(r)  n(£  ,  r)d£l  =  dr'  N(r')        -* d£    (8) 

)o     ^£,) 

where  €.  is  the  maximum  energy  at  r '  and  the  actual  number 

of  electrons  with  energy  between  £  and  c  +  d£,  namely, 

N(r')n(£',  r*)d£  has  been  introduced. 

An  exact  solution  of  this  integral  equation,  and  the  deter- 
mination of  the  distribution  function  itself,  appears  to  be  out 
of  question,  but  the  average  energy  may  be  obtained  by  the  fol- 
lowing approximate  method. 

It  can  be  seen  that  n(£! ,  r)  has  the  dimensions  of 
(energy)"1,  and  hence  can  be  written  as 

n<£,  r)  =U—Jt(e/£kV)  (9) 
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where  0   is  a  dimensionless  function  of  the  indicated  dimension- 
less  argument  and  the  dependence  on  r  is  contained  in  £^y(r). 
The  essential  approximation  consists  in  the  fulfillment  of  the 
integral  equation  only  for  those  energies  for  which  the  number 
of  electrons  is  greatest.   Since  more  than  half  of  the  electrons 
at  r  were  formed  in  the  last  mean  free  path,  and  since  the  more 
remote  the  point  of  origin  (i.e.,  the  greater  the  energy), the 
energy  distribution  will  be  monotonically  decreasing  with  the 
most  important  energy  region  at  6=0.   Thus  by  setting  r '  =  r, 

€.  »  o, 

d£     d<£-   c  V0 

dr '     dr     r 
The  average  energy  can  be  obtained  at  once  from  equation  (8) 
with  the  help  of  equations  (6)  and  (9): 

V(o)  cv0\* 

e  (r)  s  (10) 


aNm  r 


Prom  the  above  equation  we  see  that  the  average  energy,  and 
the  width  of  the  energy  distribution  which  will  be  proportional 
to  £^V'  D°th  increase  as  the  electrons  approach  the  wire,  as  was 
to  be  expected.   Prom  equation  (6)  the  average  number  of  ionizing 
collisions  per  unit  distance  follows  immediately 

K(r)  =  (0(0)  aNm  C  V^r)*  (11) 

The  final  step  in  the  calculation  is  the  consideration  of 
r0.   According  to  the  way  in  which  this  quantity  was  introduced 
in  equation  (2),  it  is  the  distance  from  the  wire  to  the  point 
at  which  the  avalanche  starts,  and  at  this  point  one  may  expect 
fluctuations  to  play  a  more  important  role.   Calculation  of  this 
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constant  in  terms  of  the  physical  constants  will  give  an  expres- 
sion for  multiplication. 

According  to  equation  (10),  introducing  a  factor  £  ,  which 
is  defined  as  the  "critical  energy",  and  is  the  average  energy 
at  the  point  where  the  ionization  starts,  can  be  written  as 

/>(0)   (C  V0)  \ 

6o" -J  (12) 

V*  Nm    ro  / 

The  critical  energy  £ Q   will  depend  on  the  nature  of  the 
gas  (ionization  potential  and  cross  section)  but  will  not  depend 
on  the  counter  voltage  or  geometrical  factors.   It  follows  from 
equation  (12)  that  r0  increases  linearly  with  the  counter  volt- 
age for  which  ionization  does  not  begin  until  the  electron 
reaches  the  wire,  i.e.,  r0  =  r^.   At  this  voltage  for  the  pro- 
portional region,  the  multiplication  factor  becomes  unity  and 
below  this  voltage  the  counter  acts  as  an  ionization  chamber. 

Denoting  the  proportional  threshold  by  Vp,  the  expression 
for  r0  can  be  written  as 

r0  =  *1  VVp  (13) 

With  equations  (13)  and  (11),  the  expression  for  multiplication 
can  be  written  as 

A  =  exp  2(aNm  C  rx  0(0)    VQ)*  [(V^)*  -  ll       (i4) 
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in  which,  as  mentioned  above,  0(0)  is  to  be  set  equal  to  unity 
and  Vp  is  to  be  determined  from  the  measurements.   The  quantity 
a  in  equation  (14)  is  the  rate  of  change  of  the  ionization  cross 
section  for  energies  just  above  the  ionization  potential. 


■'■The  quantity  0(0)  which  occurs  in  these  results  and  which 
is  essentially  the  number  of  slowest  electrons,  cannot  be  evalu- 
ated uniquely  without  knowing  the  actual  form  of  the  distribu- 
tion function.   However,  the  function  0  is  sufficiently  re- 
stricted so  that  0(0)  cannot,  in  general,  vary  over  wide  limits. 
In  addition  to  the  requirement  that  0  shall  be  a  monotonically 
decreasing  function  of  its  argument,  there  are  two  further  con- 
ditions to  be  imposed  on  this  function.   With  x  =  £/€■  AV,  the 
normalization  condition  is 

/-co 
0(x)  dx  =  1 

and  from  the  definition  of  x, 


I 


OS 


x  0(x)  dx  =  1 


Various  functional  forms  for  0  may  be  assumed  and  after  the 
first  condition  has  been  imposed,  0  may  be  calculated. 
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A  glass  vacuum  system  was  constructed  for  purification  of 
boron  trifluoride  gas  and  filling  of  the  BF3  counters.   Pressures 
of  the  order  of  10""  mm  Hg  were  achieved  and  could  be  maintained 
for  intervals  of  24  hours  by  this  system. 

Two  different  counters  (A)  and  (B)  were  constructed,  cleaned 
thoroughly,  sealed  to  the  vacuum  system,  outgassed  at  250°  C  for 
12  hours  and  at  400°  C  for  two  hours,  with  the  pressure  main- 
tained below  10"'  mm  Hg.   The  counters  were  filled  with  BF-j  gas 
at  32.0  and  29.9  cm  pressures,  respectively.   Details  of  con- 
struction of  the  counters  and  the  effect  of  impurities  in  BF3 
gas  on  the  counter  performance  are  discussed. 

Boron  trifluoride  gas  was  prepared  by  thermal  dissociation 
of  96  per  cent  B10  enriched  BF^CaFg  and  purified  by  repeated 
condensation  at  liquid  air  temperatures.   Purity  of  the  gas  was 
checked  by  measuring  its  boiling  point,  using  a  Cu-Constantan 
thermocouple.   Typical  plots  of  the  time-temperature  variation 
of  the  boiling  BF3  are  given. 

Variations  of  gas  multiplication  of  counter  (B)  with  changes 
in  the  voltage  applied  to  the  counter  were  studied.   Experimental 
values  of  multiplication  were  compared  with  their  corresponding 
theoretical  values  calculated  from  equations  derived  by  Rose 
and  Korff  (20).   Discrepancies  between  the  theoretical  and  ex- 
perimental values  at  lower  voltages  were  attributed  to  the 
asymptotic  nature  of  the  Rose  and  Korff  equation. 

Sensitivity  of  counter  (B)  for  thermal  neutrons  was  deter- 
mined by  inserting  it  in  several  positions  in  the  Kansas  State 
University  graphite  subcritical  reactor  which  has  been 
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standardized  so  that  the  neutron  flux  is  known  throughout  the 
assembly.   Dependence  of  the  counter  sensitivity  on  the  distance 
between  the  counter  and  the  source  was  studied.   Sensitivity  was 
found  to  increase  with  the  increase  in  distance  up  to  about  1.5 
diffusion  lengths  where  it  remained  constant.   This  was  attrib- 
uted to  the  fact  that  at  smaller  distances,  neutrons  are   not 
adequately  slowed  down  and  the  cross  section  of  B10  (n,  <^)   re- 
action is  correspondingly  lower.   Experimental  values  of  sensi- 
tivity were  always  much  less  than  the  theoretical  sensitivity 
calculated  by  the  formula  given  by  Price  (19).   Discrepancy  was 
attributed  to  eccentricities  in  anode  wire,  incomplete  cleaning 
of  the  counter  and/or  inaccurate  evaluation  of  the  counter  sen- 
sitive volume. 

A  number  of  plateau  curves  were  obtained  for  different 
amplifier  gain  settings.   The  best  plateau  was  obtained  at  a 
gain  of  8.   Increase  of  amplifier  gain  decreased  the  horizontal 
part  of  the  plateau. 

An  integral  bias  curve  was  obtained  applying  1475  volts  to 
the  counter.   The  relatively  broad  distribution  in  the  number  of 
counts  was  attributed  to  the  partial  ionization  loss  of  the  dis- 
integration products  which  strike  the  wall  of  the  counter. 

The  efficiency  of  the  counter  to  gamma  rays  relative  to  its 
efficiency  to  neutrons  was  determined  at  different  amplifier 
gain  settings.   This  relative  efficiency  was  found  to  increase 
with  increase  of  amplifier  gain  but  not  to  exceed  5  x  10"4. 

The  diffusion  length  of  neutrons  in  graphite  was  determined 
in  the  Kansas  State  University  graphite  subcritical  assembly 


using  counter  (B).   Based  on  one  set  of  data,  a  value  of  57.32 
cm  was  obtained.   A  better  value  could  have  been  obtained  by 
improving  the  statistics  of  counting,  by  repeating  the  experi- 
ment, and  by  allowing  for  the  flux  depression  due  to  the  presence 
of  this  relatively  large  counter.   Since  the  purpose  of  the 
experiment  was  only  to  demonstrate  the  quality  of  the  counter, 
no  further  attempts  were  made  to  measure  diffusion  length. 


